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ABSTRACT 

We report on a 151 ks net exposure Suzaku observation of the Narrow Line Seyfert 1 
galaxy Mrk 335. The 0.5-40 keV spectrum contains a broad Fe line, a strong soft excess 
below about 2 keV and a Compton hump around 20-30 keV. We find that a model 
consisting of a power law and two reflectors provides the best fit to the time-averaged 
spectrum. In this model, an ionized, heavily blurred, inner reflector produces most of 
the soft excess, while an almost neutral outer reflector (outside ^ AQ rg) produces most 
of the Fe line emission. The spectral variability of the observation is characterised by 
spectral hardening at very low count rates. In terms of our power-law -I- two-reflector 
model it seems like this hardening is mainly caused by pivoting of the power law. 
The rms spectrum of the entire observation has the curved shape commonly observed 
in AGN, although the shape is signiflcantly flatter when an interval which does not 
contain any deep dip in the lightcurve is considered. We also examine a previous 133 ks 
XMM-Newton observation of Mrk 335. We find that the XMM-Newton spectrum can 
be fitted with a similar two-reflector model as the Suzaku data and we confirm that 
the rms spectrum of the observation is flat. The flat rms spectra, as well as the high- 
energy data from the Suzaku PIN detector, disfavour an absorption origin for the soft 
excess in Mrk 335. 

Key words: galaxies: active - galaxies: Seyfert - galaxies: individual: Mrk 335 - 
X-rays: galaxies 



1 INTRODUCTION 

The X-ray spectra of type 1 Active Galactic Nuclei (AGN) 
are usually dominated by a power-law component. They also 
often exhibit a broad, skewed Fe line around 6 keV, as well 
as smooth component rising above the power-law continuum 
below about 2 keV, the so called soft excess. While it is well 
established that the broad Fe line arises due to reflection of 
the continuum power law in the inner accretion disc (where 
it is broadened by Doppler and relativistic effects) , the origin 
of the soft excess remains a question of much debate. 

The soft excess has a very similar spectral shape in 
AGN covering several decades in mass (Gierlihski & Done 
2004; Crummy et al. 2006), suggesting that its origin is in 
atomic rather than thermal processes. An obvious candidate 
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of atomic origin is the reflection spectrum associated with 
the broad Fe line. In this model, the soft excess is due to 
the blurring of soft emission lines from the inner parts of 
the accretion disc. The reflection model has been applied 
successfully in a number of AGN (e.g. Crummy et al. 2006), 
although sources that exhibit a strong soft excess but no 
broad Fe line are diflicult to explain in this model (see e.g. 
Brenneman et al. 2007). Another proposed explanation for 
the soft excess is that it is due to smeared absorption by op- 
tically thin material in the line of sight (Gierlihski & Done 
2004) . This model has been seen to provide equally good flts 
as the reflection model (e.g. Middleton et al. 2007). However, 
the model used in these flts has been very simplistic, ignoring 
emission resulting from the absorber and the acceleration of 
the matter. A more realistic absorption model does not re- 
produce the smooth shape of observed soft excesses unless 
extreme velocities are incorporated (Schurch & Done 2007). 
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Mrk 335, also known as PG 0003+199, is a Narrow Line 
Seyfert 1 (NLSl) galaxy at redshift z = 0.026. It exhibits 
both a strong soft excess and a broad Fe line. Unlike most 
Seyfert 1 galaxies, it does not show any clear signs of com- 
plex, warm absorption at low energies, making it an ideal 
target for studying the origin of the soft excess. 

The soft excess in Mrk 335 was first observed by EX- 
OSAT (Pounds et al. 1986; Turner & Pounds 1988) and was 
later confirmed by BBXRT (Turner et al. 1993). No signifi- 
cant absorption edges were detected in the BBXRT observa- 
tion but an Fe edge in Ginga data lead Turner et al. (1993) to 
suggest the presence of a variable, ionized absorber. A sub- 
sequent ASCA observation again revealed the presence of a 
strong soft excess, but showed no clear evidence of edges or 
ionized absorption at low energies (Reynolds 1997; George 
et al. 1998). The ASCA spectrum also showed evidence for 
a broad Fe line (Nandra et al. 1997), and Ballantyne et al. 

(2001) found a good fit to the entire 0.6-10 keV band with 
an ionized reflection model. BeppoSAX data (Bianchi et al. 
2001) clearly showed a broad Fe line as well as a soft excess, 
but revealed a rather small Compton reflection component. 

XMM-Newton observed Mrk 335 for the first time in 
2000. The observation was first analysed by Gondoin et al. 

(2002) , who favour a model in which the soft excess consists 
of a combination of Bremsstrahlung emission and ionized 
reflection from the accretion disc. No absorption edges apart 
from the 0.54 keV edge due to Galactic Oxygen were seen 
in the RGS spectrometer. 

The XMM-Newton observation was later re-analysed by 
Longinotti et al. (2007), who detected a narrow absorption 
feature at 5.9 keV, which, if it is identified with Fe xxvi, is 
inflowing at ~ 0.11 — 0.15 c, or, if at rest, is located close to 
the black hole and gravitationally redshifted. Both Crummy 
et al. (2006) and Middleton et al. (2007) have also fltted the 
XMM-Newton data as a part of a sample of type 1 AGN. 
Crummy et al. (2006) find that the spectrum can be well 
fitted with a reflection model, while Middleton et al. (2007) 
show that an equally good flt can be obtained if the soft 
excess is modelled as smeared absorption. 

In January 2006, XMM-Newton re-observed Mrk 335 for 
133 ks. In addition to a strong soft excess below 2 keV, this 
observation revealed a double-peaked Fe emission feature 
with peaks at 6.4 and 7 keV (O'NeiU et al. 2007). Spectral 
fltting of this feature suggested that a moderately broad rel- 
ativistic line is present but that the two peaks are due to 
narrow lines originating in more distant material. The au- 
thors further note that the reflection model inferred from 
the data above 3 keV still leaves a soft excess when extrapo- 
lated to lower energies. The rms variability spectrum of the 
observation is consistent with being constant, disfavouring 
the absorption model for the soft excess, which, assuming 
that the ionization state of the absorber is driven by the 
continuum variability, predicts enhanced variability at low 
energies (Gierlihski & Done 2006). 

Recent Swift observations of Mrk 335, carried out in 
May and June/ July 2007, caught the source in an extremely 
low state, with the flux having diminished by a factor of 
30 compared to previous observations. In this low state the 
source exhibited a very hard spectrum above 2 keV as well as 
a soft excess at low energies. The spectral changes together 
with the sudden drop in flux have been interpreted in terms 



of a partial covering model (Grupe et al. 2007) as well as in 
a reflection model (Gallo et al. 2007, in prep). 

In this paper we present the results of a 316 ks (~ 150 ks 
net exposure) Suzaku observation of Mrk 335 performed in 
June 2006. This observation was carried out before the sud- 
den drop in flux mentioned above and caught the source 
in its typical flux state. We also analyse the 133 ks XMM- 
Newton observation performed in January 2006. 

This paper is organised as follows: Section [2] describes 
the observations and the data reduction. Sections [3] and 
13] describe the spectral analysis of the Suzaku and XMM- 
Newton time-averaged spectra, and the spectral variability 
of both observations is presented in Section[5] Section [6] con- 
tains a discussion of the results, and a summary is given in 
Section [T] 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 The Suzaku observation 

Mrk 335 was observed by Suzaku between 2006 June 21-24 
for a total duration of 316 ks. Event flies from version 1.2.2.3 
of the Suzaku pipeline processing were used and spectra were 
extracted using xselect. 

The net exposure time of all four XIS detectors is 151 ks. 
For each XIS, source spectra were extracted from circular re- 
gions of 4.3 arcmin radius centred on the source (which was 
observed off-axis in the HXD nominal position). Background 
spectra were extracted from two circular regions with the 
same total area as the source region, avoiding the chip cor- 
ners with the calibration sources. Response matrices and an- 
cillary response files were generated for each XIS using xis- 
RMFGEN and xissiMARFGEN version 2006-11-26. The ARF 
generator should account for the hydrocarbon contamina- 
tion on the optical blocking filter (Ishisaki et al. 2007). 

The source count rates over the 0.7-10 keV energy range 
for the three front-illuminated (FI) detectors are 1.299 ± 
0.003 counts s-^(XISO), 1.484±0.003 counts s"^ (XIS2) and 
1.364 ±0.003 counts s"^ (XIS3). The back-illuminated (BI) 
detector (XISl) has a count rate of 2.243 ± 0.004 counts s^^ 
in the 0.3-8 keV band. The background typically contributes 
2-3 per cent of the total count rate for all four detectors. 

At the time of our observation, the bias voltage for 16 
of the 64 HXD PIN diodes had been reduced from 500 V to 
400 V, in order to suppress a rapid increase of noise events. 
The PIN event file was therefore filtered to exclude the 
diodes biased with 400 V, and the corresponding response 
file (ae_hxd_pinhxnoml23_20060814.rsp) was used. A model 
for the non- X-ray background was provided by the HXD 
team. Source and background spectra were constructed from 
identical good time intervals and the exposure time of the 
background spectrum was increased by a factor of 10 (to ac- 
count for the fact that the background model was generated 
with 10 times the actual count rate in order to minimise the 
photon noise). After deadtime correction the net exposure 
time of the PIN was 120 ks. 

The total PIN count rate over the 14-40 keV energy 
range is 0.296 ± 0.002 counts s~\ compared to 0.262 ± 
0.005 counts s~^ for the background. Since the background 
is so much higher than the source, the accuracy of the back- 
ground model is very important for our results. We therefore 
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compared the background model to the night earth spec- 
trum, which was obtained by selecting periods of earth oc- 
cultation during the observation. Since the earth is dark in 
the hard X-rays the background model should agree with 
the night earth spectrum. We find that that the background 
model overpredicts the night earth data by about 6 per cent, 
but that the spectral shape is well reproduced. We can there- 
fore simply correct the count rate of the background model 
so that it matches the night earth data when fitting the 
time-averaged spectrum. Investigations of the background 
model on shorter time-scales unfortunately show that the 
systematic uncertainties of the model are too large for tim- 
ing analysis to be carried out, and we therefore defer this 
analysis to a later time when a more accurate background 
model will be available. 

Since the background model discussed above does 
not include the contribution from the cosmic X-ray back- 
ground (CXB), a spectral model for it is included in 
all the fits. The model, which has the form 2.06 x 
10"''(S/100 keV)"^-^'^ exp(-E/41.13 keV), is based on the 
HEAO-Al spectrum and has been renormalized to the HXD 
field of view. The CXB flux over the 14-40 keV band pre- 
dicted by this model represents about 25 per cent of the 
source flux. 

2.2 The XMM-Newton observation 

XMM-Newton observed Mrk 335 between 2006 January 
03-05 for a total duration of 133 ks (observation ID 
0306870101). The data were reduced using the XMM- 
Newton Science Analysis System version (7.0.0). We use only 
the data from the EPIC pn camera, which was operated in 
the Small Window mode during the observation. Because of 
the Small Window mode the data are not affected by photon 
pile-up. The source spectrum was extracted from a circular 
region of radius 34 arcsec centred on the source, and the 
background spectrum was extracted from two circular re- 
gions with a total area three times that of the source region. 
Inspection of the background showed that strong flaring was 
present at the beginning and at the end of the observation. 
Excluding these time intervals leaves a total good exposure 
time of 115 ks. 



3 THE TIME-AVERAGED SUZAKU 
SPECTRUM 

In order to investigate the agreement between the four XIS 
detectors we initially fitted them separately with a simple 
power-law model. The three FI XIS showed very similar 
spectra over the 0.7-10 keV energy range, with the photon 
indices agreeing to within 0.01. Larger differences are seen 
below 0.7 keV and the agreement becomes very poor below 
about 0.6 keV (see e.g. Fig. [5}. Since we find that including 
the data between 0.6 and 0.7 keV only makes the quality 
of the fits worse without significantly changing any of the 
fit parameters, we choose 0.7 keV as our lower energy limit. 
(We will use the BI XIS to constrain our spectral models at 
lower energies.) When the 0.7-10 keV FI XIS spectra were 
fitted simultaneously with the photon index tied, the rela- 
tive normalizations were found to be XIS0/XIS2 = 0.93 and 
XIS3/XIS2 = 0.99. These values were also obtained when 



fitting more complicated models, and in all subsequent fits 
we will fit the FI XIS simultaneously with gamma tied. The 
BI XIS (XISl) was found to have a somewhat steeper spec- 
trum than the other XIS (AF « 0.1 in the same energy 
range). As XISl also has a higher effective area at softer 
energies (we use the 0.3-8 keV energy range) we choose to 
fit it separately and use it mainly to constrain our models 
at very low energies. Due to calibration uncertainties near 
the instrumental Si K edge at 1.84 keV we always ignore 
1.8-1.9 keV in the FI XIS and 1.7-1.9 keV in the BI XIS. 

The PIN spectrum is unfortunately associated with 
uncertainties regarding the background modelling and the 
cross normalization with respect to the XIS detectors. Al- 
though this means that we cannot use the PIN data to ac- 
curately constrain parameters of different spectral models, 
we can still use it to distinguish between models that fit the 
data at lower energies. 

Below we will first consider the 2-10 keV spectrum and 
then include the low-energy data down to 0.3 keV. As a last 
step we will consider the 14-40 keV PIN spectrum. Through- 
out this paper, errors are quoted at the 90 per cent confi- 
dence level and energies of spectral features are quoted for 
the rest frame of the source. 



3.1 The 2-10 keV spectrum 

The most prominent feature in the 2-10 keV energy range 
is the broad Fe line which extends from about 5 to 7 keV. 
The line profile is shown in Fig. [1] as a ratio to a power law 
modified by Galactic absorption (Nh = 3.99 x 10^" cm~^, 
included in all fits from hereon), fitted over the 2-4.5 and 
7.5-10 keV energy ranges. The power law has a photon index 
of r = 2.09 ± 0.02. 

In order to find a phenomenological description for the 
data over the whole 2-10 keV band we start by adding a 
Gaussian line to the power-law model. This simple model 
provides a good fit to the data with — 4441 for 4353 
degrees of freedom (d.o.f.). The Gaussian line is found to 
have E = 6.43 ± 0.04 keV, cr = 0.40t^ |;^ keV and equiv- 
alent width EW = 285 ± 34 eV. The photon index of the 
power law is still F = 2.09 ± 0.01. When adding a narrow 
line around 6.4 keV to the model (to check for a narrow 
component of the Fe line from distant matter) the fit im- 
proves by Ax^ = 19 for 2 de grees of freedom. The nar- 
row line has energy E = 6.36 ± 0.03 keVQ and equivalent 
width EW = 29 ± 14 eV. With the narrow line included, 
the broad line parameters change to E = 6.45 ± 0.05 keV, 
o = 0.47t° J° keV and EW = 250t^° eV. 

No other narrow lines are detected in the spectrum. 
In particular, we do not detect the 5.9 keV absorption line 
seen in the 2000 XMM-Newton observation (Longinotti et 
al. 2007) or the 7 keV emission line seen in the 2006 XMM- 
Newton observation (O'Neill et al. 2007). If the narrow line 
around 7 keV originates in distant material, as suggested 



^ After completion of this work a new processing version of these 
data was released (V. 2. 0.6. 13). With this version the energy of 
the narrow line is found to be E = 6.391q °3 keV, i.e completely 
consistent with the 6.4 keV Fe line. No other model parameters 
changed significantly with the new processing version. 
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Figure 1. The Suzaku FI XIS spectrum shown as a ratio to a 
power law fitted over the 2-4.5 and 7.5—10 keV energy ranges. 
Data points in black, red and green correspond to XIS 0, 2 and 
3. 

by O'Neill et al. (2007), it seems strange that is has dis- 
appeared in the 6 months between the XMM-Newton and 
Suzaku observations. However, adding a line around 7 keV 
to our model for the Suzaku data does not improve the qual- 
ity of the fit, and we find that the line tends to move to 
energies much lower than 7 keV. If we freeze the line en- 
ergy at 6.97 keV (the energy of Fe xxvi Lya) we find a 
line flux of O.SSlgg x 10~®photons cm~^ s~^, corresponding 
to an equivalent width of EW = 6lg^ eV. Using the same 
model for the 2006 XMM-Newton data, we find a signifi- 
cantly higher line fiux of 4.96tl'82 ^ 10~^photons cm~^ s~^. 
It thus seems clear that the emission line around 7 keV 
has disappeared or at least weakened significantly in the 
Suzaku observation. 



We now move on to try and find a more realistic model 
for the data. The broad, skewed Fe line seen in Fig. [1] arises 
as the incident power law is reflected in the accretion disc, 
and should therefore have a reflection continuum associated 
with it. To model the reflection from an accretion disc we 
use the constant-density model reflion by Ross & Fabian 
(2005). The parameters of the model are the Fe abundance, 
the ionization parameter ^, the photon index of the inci- 
dent power law and the normalization. In order to account 
for relativistic effects in the vicinity of the black hole, this 
model is convolved with the relativistic blurring kernel kd- 
BLUR, which is derived from the code by Laor (1991). The 
relativistic blurring parameters are the inner and outer radii 
of the disc (rin and rout), the inclination i and the emissivity 
index q (the emissivity follows the form e oc r~', where r is 
the radius of emission). 

In addition to the power law and the blurred reflection 
component, we also include the narrow Fe Ka line in our 
model. In all flts we fix the Fe abundance at the solar value, 
the emissivity index at g = 3 and the outer radius of the 
disc at Tout ~ 400 rg (where rg — GM/c? is the gravita- 
tional radius). We find a best fit of xVd-o.f. = 4416/4353 
with r = 2.19 ± 0.01, i = 30+^ erg cm s"^ (the lowest al- 
lowed value), i — 54° and rin — 40I20 ''s- parameters 
of the narrow line are found to be E = 6.36 ± 0.02 keV and 
EW = 60^1c eV. In order to confirm these results we also 
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Figure 2. The Suzaku XIS spectrum shown as a ratio to the 2— 
10 keV power-law -|- Gaussian model, clearly revealing the soft 
excess. Data from the FI XIS are shown in black, red and green, 
and data from the BI XIS (XISl) are shown in blue. The agree- 
ment between the detectors clearly becomes poor at low energies. 
The structure seen in the 0.5-0.7 keV range in XISl is discussed 
in the text. 

fitted a model in which we replaced the blurring kernel kd- 
BLUR with the kbrrconv model by Brenneman & Reynolds 
(2006). An important difference between these two models 
is that the kerrconv model has the spin of the black hole 
as a free parameter. The kerrconv model gives best-fitting 
parameters that are consistent with those reported above, 
and gives a best-fitting spin value of 0, as expected from the 
relatively large value of rin. 

It should be noted that the best fit presented above is 
not necessarily unique. Both a highly ionized disc and a high 
inclination makes the blue peak of the Fe line bluer, and 
these parameters are hence somewhat degenerate. In fact, 
with an ionisation parameter of ^ = 2000 erg cm s"'^, we 
find a fit of comparable quality to the one above (x^/d.o.f . — 
4424/4353) with i = 19° IJg, F = 2.06 ± 0.01 and similar 
parameters for the narrow line. We will hopefully be able to 
break this degeneracy by including the low-energy data. 

3.2 The soft excess 



We now include the data from the three FI XIS down to 
0.7 keV (this lower energy limit is discussed at the begin- 
ning of the section), ignoring the 1.8-1.9 keV band which 
is affected by the instrumental Si edge. When extrapolating 
the best-fitting 2-10 keV models, a clear, rather smooth soft 
excess is visible below 2 keV, as shown in Fig. [21 The ratio 
of the data to the model at 0.7 keV is 1.9 for the simple 
power-law model, and 1.7 and 1.5 for the highly ionized and 
almost neutral reflection models, respectively. 

A good fit to the soft excess can be found by simply 
adding a black-body component to the 2-10 keV models. 
A black-body of temperature 0.13 keV together with our 
reflection model with ^ = 30 erg cm s~^ for instance gives 
X^/d.o.f. — 5579/5340. Although this solution provides a 
good parametrisation of the data it is very unlikely to be 
realistic. The temperature is simply too hot for an accretion 
disc around a supermassive black hole, and black-body tem- 
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Figure 3. Best-fitting two-refiector model for the time-averaged 
0.7—10 keV Suzaku spectrum. The solid, black line shows the total 
model and the dashed, red line shows the primary power law. The 
inner and outer reflectors are represented by the blue, dotted line 
and the green, dashed-dotted line, respectively. The model also 
includes a narrow Fe Ka line. 
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Figure 4. The best-fitting blurred absorption model for the 
time-averaged 0.7—10 keV Suzaku spectrum. The solid, black line 
shows the total model and the dashed, red line shows the absorbed 
power law. The model also includes distant refiection (green, 
dashed-dotted line) and a narrow Fe K« line. 



peratures of 0.1-0.2 keV have been found to fit soft excesses 
in AGN with a large range of masses and accretion rates 
(Gierliiiski & Done 2004; Crummy et al. 2006), in disagree- 
ment with the prediction that the disc temperature should 
scale as T oc M-^/'^{M / MEddY'^ (Shakura & Syunyaev 
1973). 

A natural explanation for the constant temperature of 
the soft excess is that it is due to some atomic process, and 
reflection is thus an obvious candidate. Since no emission 
lines are seen in the soft excess, this reflection has to be 
heavily blurred, and therefore originate in the very inner 
parts of the accretion disc. This poses a problem for our 
2-10 keV reflection models. The inner radius of emission at 
40 Tg, as suggested by the Fe line shape, simply does not pro- 
vide enough blurring to reproduce the smooth shape of the 
soft excess. A disc truncated at ~ 40 rg is also at odds with 
Mrk 335 being a high accretion rate source (M/MEdd = 0.79, 
based on Woo & Urry 2002) . It is thus possible that the disc 



extends further in but is unable to produce fluorescent Fe 
emission in the inner parts (e.g. because of its ionization 
state). 

It seems plausible that a model with two reflectors can 
reproduce the data; an inner, heavily blurred reflector to 
produce the soft excess, and an outer reflector to produce 
the sharp features of the Fe line. We therefore construct a 
model with an inner reflector extending from from 1.24 to 
40 Tg and an outer reflector extending from 40 to 400 r-g. If 
we fix the inclination at i = 54° (as found for our best fit to 
the 2-10 keV data) we find a fit with xVd-o.f. = 5592/5335. 
This model has a photon index of F = 2.19 ± 0.01 and inner 
and outer ionization parameters of ^in = 310lg'' erg cm s~^ 
and ^out = 321^2 ^rg cm s~^. The narrow line is found to 
have energy E = 6.36 ±0.02 keV and an equivalent width of 
EW = 67^ J4 eV. The reflection fraction of the two refiectors 
is 2.8, with each of the reflectors contributing roughly the 
same amount. The model is plotted in Fig. |3] 

It is clear that this three-component model leaves much 
freedom for many of the parameters and that the fit pre- 
sented above is not necessarily unique. The fact that the 
inner disc needs to be more highly ionized than the outer 
one however seems robust. We were for example not able to 
find a satisfactory fit with a highly ionized outer reflector 
(together with a lower inclination) as suggested by one of 
our flts to the 2-10 keV data. 

Another possible way of producing the soft excess is 
by relativistically blurred absorption of a steep continuum. 
Just as for the case of reflection, this model explains both 
the constant temperature and the smooth shape of the soft 
excess. In order to investigate this possibility we use the 
swiNDi model of Gierlinski & Done (2006). This model has 
three free parameters, the velocity smearing a (assumed to 
be Gaussian), together with the ionization parameter and 
the column density of the absorbing material (^a and Nh re- 
spectively). To model the continuum we use a power law and 
a blurred reflection component, as the latter is still required 
to explain the broad Fe line. For the reflection component 
we use the same parameters as those found for the outer 
reflector in the two-reflector model, leaving only the nor- 
malization as a free parameter. We find a fit of comparable 
quality as the two-reflector model (x^/d.o.f. — 5581/5334) 
with = 499tg erg cm s-\ Nh = 7.4 ± 0.9 x 10^^ cm'^, 
a = 0.50-0.02 c (highest allowed value) and F = 2.47to;o2- 
The model is shown in Fig. |4l 

The fact that the soft excess can be equally well fitted 
with smeared absorption as with blurred refiection from an 
accretion disc is commonly observed in AGN (e.g. Middleton 
et al. 2007). However, the velocity smearing of 0.5 c that we 
find in our best fit with the absorption model is unreason- 
ably high. It should also be noted that the swindi absorp- 
tion model is rather simplistic in its treatment of the velocity 
field (assumed to be Gaussian) and in that it does not in- 
clude any emission from the absorbing material. It has been 
shown that when these issues are addressed, the model spec- 
tra include sharp features which contrast strongly with the 
observed smooth shapes of soft excesses (Schurch & Done 
2007). We therefore favour the two-reflector model presented 
above as an explanation for the soft excess in Mrk 335. 
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3.2.1 The spectrum from the back-illuminated XISl 
detector 

For XISl, we use the 0.3-8 keV energy band, ignoring the 
1.7-1.9 keV range due to the effects of instrumental Si. Since 
the calibration of XISl is still somewhat uncertain we start 
by comparing its low-energy spectrum with that from the 
long XMM-Newton observation. When taking the ratio of 
the two data sets, clear, sharp differences can be seen at low 
energies. The structure seen in the 0.5-0.7 keV range in Fig. 
[2] is for instance not present in the XMM-Newton data. As 
it is quite likely that at least some of these differences are 
due to calibration uncertainties in XISl, we will only very 
briefly discuss the XISl spectrum. 

The two-reflector model provides a good fit to the XISl 
data down to 0.5 keV (x^/d.o.f. = 1614/1605), with none of 
the parameters changing significantly from the values pre- 
sented above. Extrapolation of the model down to 0.3 keV 
however reveals a soft excess as well as the very sharp edge 
around 0.5 keV. In order to fit these features with the two- 
reflector model we need to add an edge at 0.49 keV, and 
allow for a steeper emissivity proflle (g — 6.5) as well as a 
steeper power law (F = 2.29). If the 0.49 keV edge is intrin- 
sic to the source it can be identifles with C vi, although it 
seems more likely that it is an effect of inaccurate modelling 
of the contamination in XISl. 



3.3 The 14 - 40 keV PIN spectrum 

In addition to a soft excess and a broad Fe line, reflection 
from an accretion disc should produce a Compton hump 
around 20-30 keV, as seen in Fig. [3] Data at these energies, 
provided by the Suzaku HXD PIN detector, are therefore 
crucial for testing the models derived from the XIS data. 
Speciflcally, the two-reflector model predicts about 1.4 times 
more flux in the 14-40 keV range than the smeared absorp- 
tion model. 

Mrk 335 is detected in the PIN up to about 40 keV and 
we stress that this detection is robust against uncertainties 
in the background model. However, since the background 
dominates the spectrum, the accuracy of the background 
modelling is very important for our results. As mentioned in 
section [2l comparison of the instrumental background model 
with night earth spectra shows that the background model is 
overestimated by about 6 per cent. We therefore correct for 
this before fitting the spectrum. We also include a model for 
the CXB as described in the same section. The 14-40 keV 
flux from the CXB model represents about 25 per cent of 
the source flux. The cross normalization of the PIN with 
respect to the XIS detectors has been reported to be 1.16 
for the HXD nominal position and vl.2.2.3 data (Ishida et 
al. 2006), and we use this value in the following analysis. 

Fig. [5] shows the entire 0.7-40 keV spectrum as a ratio 
to the simple phenomenological 2-10 keV model discussed 
in section [3Tl1 In addition to the soft excess we see clear ex- 
cess emission in the PIN detector. In order to check if this 
emission can be reproduced by the Compton hump of our 
0.7-10 keV two-reflector model we extrapolate this model 
to 40 keV. This gives x^/d.o.f. = 5690/5410 for the entire 
0.7-40 keV band. Fitting of the model with the PIN data 
included does not improve the quality of the fit or change 
any of the parameter values. The fit is shown in Fig. [6] and 
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Figure 5. The broad-band Suzaku spectrum of Mrk 335 shown 
as a ratio to the simple 2—10 keV power law + Gaussian model. 
Excess emission is clearly visible both at low and high energies. 
Data from the three FI XIS are shown in black, red and green. 
The PIN data are shown in blue. 
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Figure 6. The broad-band Suzaku spectrum of Mrk 335. The 
upper panel shows the data together with the two-reflector model 
shown in Fig. |3] The lower panel shows the residuals of the fit. 
Data from the three FI XIS are shown in black, red and green. 
The PIN data are shown in blue. 



we see that the model slightly underpredicts the PIN data. 
Given that the source is not very bright in the hard X-rays, 
this slight discrepancy is likely to be due to uncertainties 
in the background models (both non-X-ray background and 
the CXB) and/or the cross- normalization of the PIN with 
respect to the XIS detectors. If we instead extrapolate the 
blurred absorption model to include the PIN data, we get 
X^/d.o.f. = 5734/5410, i.e Ax^ = 43 worse than for the 
two-reflector model. Fitting of the absorption model with 
the PIN data included results in a slightly higher contribu- 
tion from the reflector, which improves the quality of the fit 
by Ax^ = 4. The high -energy data hence favour the two- 
reflector model, although both models are acceptable given 
the uncertainties mentioned above. 
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Figure 7. The XMM-Newton pn spectrum shown as a ratio to 
a power law fitted over the 3-4.5 and 7.5—10 keV energy ranges. 
The Fe Une profile is clearly double peaked, but apart from the 
blue peak it is very similar to the profile shown in Fig. [l] for the 
Suzaku data. 



4 THE TIME-AVERAGED XMM-NEWTON 
SPECTRUM 



Mrk 335 was observed for 133 ks by XMM-Newton 
about six months before the Suzaku observation, and a com- 
parison between these two observations is likely to provide 
interesting information about the source. Because XMM- 
Newton is more rehably cahbrated at low energies, we also 
expect that the XMM-Newton data will allow us to put bet- 
ter constraints on the two-reflector model at low energies. 

We start, however, by considering the 2-10 keV time- 
averaged spectrum. The most notable difference compared 
to the Suzaku data in this energy range is the very clearly 
double-peaked Fe line profile, shown in Fig. [7] Detailed fit- 
ting of the line profile has been presented by O'Neill et al. 
(2007), who find that a moderately broad relativistic line is 
present but that the two peaks at 6.4 and 7 keV are most 
Ukely due to narrow lines, identified with Fe Ka and Fe xxvi 
Lya, respectively. The two lines are suggested to originate 
in the molecular torus of AGN unification models, and in 
highly ionized gas filling this torus, respectively. Here we 
simply note that the Fe line profile can be fitted with the 
same model as the Suzaku data if a narrow line is added 
around 7 keV, in agreement with the findings of O'Neill et 
al. (2007). 

With the addition of the narrow hue around 7 keV, the 
entire 2-10 keV XMM-Newton spectrum can be very well 
fitted with the models presented for the Suzaku data. As 
an example, our model with reflection arising from outside 
40 rg in a disc with ^ = 30 erg cm s~^, gives x^/d.o.f. = 
1324/1353, with the photon index and the relative nor- 
malization of the different components changing very little 
from the Suzaku values. For the two narrow lines we find 
energies of 6.411°;°^ keV and 7.05115:2^ keV, with equiv- 
alent widths of 43 ± 13 eV and 32 ± 14 eV, respectively. 
The 2-10 keV flux determined from the model is 1.78 x 
10~^^ erg cm"'^ s~^, compared to 1.43 x 10~^^ erg cm~^ s~^ 
for the Suzaku data. The flux in the narrow Fe Ka is consis- 
tent with not having changed between the observations, with 
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Figure 8. The XIS3 lightcurve in the 0.5-12 keV band in 300 s 
bins. We also show the high and low fiux states, which were de- 
fined so that they have roughly the same number of counts. 



7-4ltl;53 X 10~® photons cm s obtained for Suzaku and 
6.58ii;99 X lO"'^ photons cm'^ s'^ for XMM-Newton. 

At low energies, larger differences between the two data 
sets can be seen. In terms of the two- reflector model, much of 
this difference can be explained by increasing the ionization 
parameter of the inner reflector from about 300 erg cm s~^ 
to about 1200 erg cm s^^ (no solution with the same ioniza- 
tion parameter for the inner reflector in both data sets was 
found). However, in order to get a good fit all the way down 
to 0.4 keV we also need to include an edge at 0.51 keV with 
r — 0.33 (as for XISl, possibly consistent with the 0.49 keV 
C VI edge) . If we let all the parameters of the model vary we 
find a best-fitting 0.4-10 keV two-reflector (-1- edge) model 
with x^/d.o.f. = 1749/1670. Apart from the higher ioniza- 
tion parameter for the inner reflector, the only parameters in 
this model that differ significantly from the Suzaku FI XIS 
fits are the photon index and the emissivity index. These are 
both steeper; P = 2.28 ± 0.01 and q = 5.9lQg, respectively 
(we froze g at 3 for the FI XIS data). The reflection fraction 
and relative contribution of the two reflectors in this model 
are very close to what we found for the Suzaku data. 



5 SPECTRAL VARIABILITY 



In Fig. |8l we show the 0.5-12 keV background- 
subtracted XIS3 hght curve from the Suzaku observation of 
Mrk 335. The figure also defines the high and low states 
which we will use to study the spectral variability. The 
source varies by a factor of about 2.5 during the observa- 
tion and the lightcurve notably contains two big fiares. 

Below, we will investigate the spectral variability of 
Mrk 335 during the Suzaku observation by considering 
hardness-ratios, the high and low states, and rms spectra. 
Because of the high uncertainties in the PIN background on 
short time-scales we limit this analysis to the XIS data. At 
the end of the section we will compare our results with the 
variability properties of the XMM-Newton data. 
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Figure 9. The mean XIS hardness-ratio as a function of the mean 
count rate. The hardness-ratio is defined as H — S/H -|- S, where 
H is the count rate in the 3-10 keV band and S is the count rate 
in the 1—2 keV band. The source clearly hardens significantly at 
low count rates. 
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Figure 10. The mean XIS hardness-ratio (black, solid line) as 
a function of time together with the re-scaled, mean XIS 0.5- 
12 keV lightcurve (red, dashed line). The source hardens signif- 
icantly during the dips in the lightcurve but softens very little 
during the two big flares. The vertical dashed lines divide the 
observation into three time intervals of equal length. We will con- 
sider the rms spectra from these three intervals in section [53] 



5.1 Hardness-ratios 



As a first, model-independent, way of characterising 
the spectral variability of the source, we calculate hardness- 
ratios, HR = H — S/H -I- S, where we take H to be the count 
rate in the 3-10 keV band and S to be the count rate in the 
1-2 keV band. The mean XIS hardness-ratio, calculated in 
orbital length bins of 5760 s, is plotted as a function of the 
mean H -I- S count rate in Fig. [9] At low count rates there 
is a clear anti-correlation between the hardness-ratio and 
the count rate, but this anti-correlation flattens out as the 
count rate increases. This curious behaviour is even more ev- 
ident in Fig. 1101 which shows the mean XIS hardness-ratio 
as a function of time, together with the re-scaled, mean XIS 
0.5-12 keV lightcurve. We see that the source hardens sig- 



nificantly during the dips in the lightcurve but that it does 
not soften by a corresponding amount during the two big 
flares. 

In Fig. 1101 we also divide the observation into three time 
intervals of equal length. The hardness ratio in the middle 
part varies less than it does in the first and the third parts 
(mainly because this part of the lightcurve does not contain 
any large dip), and in section [5]3] we will see how this leads 
to different shapes of the rms spectra for the three intervals. 

5.2 High and Low states 

In the previous section we have seen that there is an anti- 
correlation between the spectral hardness and the flux of 
the source. In order to investigate what spectral compo- 
nent (s) is(are) responsible for this behaviour we now con- 
sider flux-selected spectra. We select low-, intermediate- and 
high-state FI XIS spectra as defined in Fig. [8] With these 
selection criteria, the high and low states have about the 
same number of counts, corresponding to one fourth of the 
total number of counts each. We fit both a phenomenologi- 
cal model and the two-refiector model (described in section 
I3.2|l to each of the three states over the 0.7-10 keV energy 
band. 

The phenomenological model consists of a power law, 
a black-body to model the soft excess, and a broad and a 
narrow Gaussian line to model the broad and narrow com- 
ponents of the Fe line. We let all parameters vary apart 
from the widths of the Gaussian lines, which we freeze at 
fbroad = 0.47 keV (as found for the time-averaged fits in 
section 13. ip and cTnarrow = 1 eV. The results from the fits 
are presented in Table [1] As expected, we see clear evidence 
of spectral pivoting, with the photon index of the power 
law increasing from 2.07 in the low state, to 2.12 in the in- 
termediate state, and 2.16 in the high state. (This linear 
increase of F is simply an effect of how the flux states were 
selected. When selecting a lower low state, we do indeed see 
a much harder power-law, as expected from Figs. iQland llOH 
The temperature of the black-body is found to be roughly 
constant at 0.14 keV in all three states, while its normaliza- 
tion increases with flux. Since the luminosity of a black-body 
should scale with its temperature, this confirms that the soft 
excess is not really a black-body. For the broad component 
of the Fe line, we find that the equivalent width stays nearly 
constant, showing that this component varies together with 
the continuum. The parameters of the narrow line are not 
very well constrained, but we note that the flux of the line 
is consistent with being constant, as would be expected if 
the line originates in distant material. 

We next fit the three fiux states with the two-reflector 
model from section 13.21 We freeze all parameters apart from 
F, the energy of the narrow line, and the four normaliza- 
tions (power law, two reflectors and narrow line) at the time- 
averaged values found in section [321 The results from these 
fits are also presented in Table [1] As for the phenomeno- 
logical model, the power law steepens with fiux, with F in- 
creasing from 2.14 in the low state to 2.26 in the high state. 
Any attempts to fit the data with F frozen at the time- 
averaged value of 2.19 give fits which are significantly worse 
than those presented in the Table. 

As a measure of the relative contribution of the power 
law and the two refiectors in the different fiux states, Table[T] 
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Parameter 



Low state Intermediate state High state 





Power law + black-body 


+ broad line - 


- narrow line 




r 

PL Norm 
BB Temp 
BB Norm 

-^broad 
EWbroad 
En arrow 
FluXnarrow 

xVd.o.f. 


xlO-3 
eV 

xlO-4 
keV 
eV 
keV 

xlO~^ph cm~2 g-l 


2.07 ±0.02 
4 74+8-26 

^•'*-o.io 
134+5 
1 41+0-" 
6.54«:lt 

2911:11 

6.36 ±0.04 

c 40 + 3.87 
"^•^•^-2.37 

2818/2831 


2.12 ± 0.02 
6.42 ± 0.10 
135±4 

^- '"-0.09 

6 41+0-09 
"•^^-0.08 

2941^0 

6.39 ± 0.05 
3 19+3-10 

3870/3879 


2.16 ±0.02 

8.59 ±0.20 

142+3 

2 09+0-15 

6 42+0-13 
°-^^-0.12 

288+7g 
6.3810-- 

2693/2659 


Power law + two reflectors + narrow line 


r 

PL/Tot fluxi 

Refinner/Tot flux^ 

Refoutcr/Tot fluxi 

Enarrow 
FluXnarrow 

xVd.o.f. 


(0.7-10 keV) 
(0.7-10 keV) 
(0.7-10 keV) 
keV 

xl0~''ph cm~2 


2.14±0.02 
0.72 
0.22 
0.06 

6.37 ±0.03 
7 94+2-66 

2800/2833 


2.20±0.02 

0.74 

0.19 

0.07 

6.38 ± 0.03 
09+2.30 

"•^^-2.13 

3826/3881 


2.26±0.02 
0.75 
0.18 
0.07 

6 37+0-07 
"-'-''-0.06 
c-n+4.07 
'5.0J_4 00 

2665/2661 



Table 1. Results from fits to low-, intermediate- and high-state FI XIS spectra over the 0.7-10 keV band. The upper values are for a 
simple phenomenological model and the lower values are for the two-reflector model described in section [3.21 l Error bars are not given 
as it is currently not possible to calculate errors on fluxes of speciflc model components in XSPEC. 



gives the ratio of the flux in each of these components to 
the total 0.7-10 keV flux. We find that the outer reflector 
contributes by the same amount (6-7 per cent) in all the flux 
states, while the inner reflector contributes slightly more to 
the low state (22 per cent in the low state compared to 18 
per cent in the high state), indicating that it is somewhat 
less variable. Although this effect is quite small, we note that 
a less variable inner reflector would be expected in a picture 
where light bending is important close to the central black 
hole (e.g. Miniutti & Fabian 2004). 

Since the spectral hardening is most prominent in the 
two deep dips of the lightcurve we also extracted spectra 
in these intervals. Specifically, we selected the 12 ks (~6 ks 
exposure) with the lowest count rate in each of the dips. Fit- 
ting of these spectra are consistent with the trends discussed 
above, although the short exposure times make it difficult 
to unambiguously disentangle the contribution from the dif- 
ferent spectral components. 

From the fits to both of the models presented above it 
seems clear that the main driver behind the observed anti- 
correlation between the hardness-ratio and the fiux is a piv- 
oting power law. Both the broad Fe line and the soft excess 
appear to largely vary together with the continuum, while 
the narrow component of the Fe line is consistent with being 
constant. 



5.3 Rms spectra 



Another way of gaining insight into the variability prop- 
erties of Mrk 335 is the rms spectrum, which, for a given 
time-scale, shows the fractional variability as a function of 
energy. The techniques for calculating rms spectra are de- 
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Figure 11. Rms spectra of Mrk 335 on a 6 ks time-scale. The rms 
spectrum of the entire observation is shown in black, while rms 
spectra in red, green and blue were calculated for the first, second 
and third parts of the observation (as indicated in Fig. IIOI I . The 
middle part of the observation clearly gives a much flatter rms, 
consistent with the small variation in hardness-ratio seen in Fig. 
1 101 during this time interval. 



scribed in e.g. Edelson et al. (2002) and Vaughan ct al. 
(2003). 

Fig. Illl shows rms spectra of Mrk 335, calculated on the 
orbital time-scale of 5760 s. In order to investigate whether 
the rms spectrum changes with time, we have calculated rms 
spectra from the first, second and third parts of the observa- 
tion (the three intervals are defined in Fig. I10|l in addition 
to the rms spectrum of the entire observation. The errors 
from the Poisson noise were calculated following Vaughan 
et al. (2003). 

The rms spectrum of the entire observation (black 
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crosses in Fig. Ilip peaks around 1 keV and then decreases 
with energy, a trend which is commonly observed in AGN. 
This behaviour is often explained either in terms of a piv- 
oting power law, or in in a model consisting of a constant 
reflection component and a power law that varies in nor- 
malization but not in shape. It can also be explained in 
an absorption model, assuming that the ionization param- 
eter of the absorber is driven by the continuum variability 
(Gierlihski & Done 2006). In the case of Mrk 335, the fits to 
the different flux states in the previous section suggest that 
the shape of its rms spectrum is mainly due to a pivoting 
power law, but that a somewhat less variable inner reflector 
may also contribute. 

The rms spectra from the three parts of the observation 
show a change in shape as well as in normalization, with 
the middle part being significantly fiatter than the first and 
the last parts, which both have roughly the same shape as 
the rms of the whole observation. The flatter rms of the 
middle part is consistent with Fig. 1101 which shows that 
the hardness-ratio is less variable during this period (the 
hardness only changes significantly during the large dips in 
the lightcurve, none of which are present during the middle 
part). The change of the shape of the rms spectrum is hence 
simply due to the position and the length of the time interval 
that is being probed, and does not imply that the source has 
entered a new mode of variability or changed its spectrum 
(the time-averaged spectra of the three time intervals are 
essentially identical). We also note that the fiat shape of 
the rms is difficult to explain within the absorption model, 
which, assuming that the ionization state of the absorber 
responds to the variable continuum, predicts an enhanced 
variability in the 0.8-2 keV range (see also O'Neill et al. 
2007). 

In order to confirm that the two-reflector model is able 
to accurately describe the variability, we also produced syn- 
thetic rms spectra for the three intervals and the entire ob- 
servation. This was done by fltting the model to the individ- 
ual 6 ks spectra and then using the best-fitting models to 
extract fluxes in the different energy bands. We flnd a good 
match between the observed and synthetic rms spectra as 
long as we let the photon index and the normalizations of the 
three components be free parameters in the fits, as expected 
from the results of the previous section. Although this result 
is encouraging, it should be noted that this simply means 
that our two-refiector model provides a good fit to the entire 
energy range in all the individual time intervals. 



5.4 Comparison with the XMM-Newton 
observation 



In section |4l we have seen that the long XMM- 
Newton observation of Mrk 335 catches the source in a very 
similar flux state as the Suzaku observation does six months 
later. The time-averaged spectra of the two observations are 
very similar above 2 keV (apart from an emission line at 
7 keV only present in the XMM-Newton data), and we find 
that similar two-reflector models flt both data sets also at 
lower energies. In this section we will extend the compar- 
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Figure 12. Same as Fig. [10] but for the XMM-Newton obser- 
vation. Note that both the lightcurve (red, dashed line) and the 
hardness-ratio (black, solid line) are much less variable than in 
the Suzaku observation. 
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Figure 13. Rms spectra of the entire XMM-Newton observation 
on a 1 ks (black) and a 6 ks (red) time-scale. The 1 ks rms is 
consistent with being constant at 0.12. The 6 ks rms decreases 
slightly with energy but it is much flatter than the rms of the 
whole Suzaku observation. The inset shows the 5-10 keV range 
of the 6 ks rms with a much finer energy binning. There is a clear 
dip in the variability at 7 keV, confirming that the 7 keV peak 
seen in the spectrum originates in distant material. 



ison between the two observations by considering some of 
the variability properties of the XMM-Newton data. 

Fig. [12] shows the hardness-ratio together with the re- 
scaled lightcurve of the XMM-Newton observation, calcu- 
lated in the same time bins as for the Suzaku observation (cf. 
Fig. nop. We note that both the count rate and the hardness- 
ratio are much less variable than in the Suzaku observation, 
and that they do not show a correlation. 

From the almost constant hardness-ratio we would ex- 
pect the rms spectrum of the XMM-Newton observation to 
be fairly flat. This is exactly what we see in Fig. 1131 which 
shows the rms spectrum on both a 1 ks and a 6 ks time- 
scale (the Suzaku rms spectra were calculated on a 6 ks 
time-scale). The 1 ks rms spectrum is consistent with be- 
ing constant at 0.12, in excellent agreement with O'Neill 
et al. (2007). The 6 ks rms decreases slightly with energy, 
but is still much flatter than the rms spectrum of the en- 
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tire Suzaku observation. This difference wouid naiveiy sug- 
gest that the source has changed its behaviour dramaticaiiy 
in the six months between the two observations. However, 
we have seen that the Suzaku rms changes in shape during 
the observation, and that the curved shape arises when the 
Ughtcurve contains a iarge dip, in which the source hard- 
ens significantiy. It therefore seems plausible that the rms 
spectrum is flat during the XMM-Newton observation sim- 
ply because the count rate never drops to a low enough level. 
This interpretation is further supported by the fact that the 
flux is slightly higher in the XMM-Newton observation and 
that it is only a third as long as the Suzaku observation. 

Fig. [13] also shows the 5-10 keV range of the 6 ks rms 
with a much finer energy binning, revealing a clear dip in 
the variability at 7 keV. This supports our interpretation 
that the 7 keV peak seen in the spectrum is a narrow line 
originating in distant material 



6 DISCUSSION 

Mrk 335 exhibits a strong soft excess below about 2 keV 
and does not show any signs of complex, warm absorption 
at low energies. This makes it an ideal target for studying 
the origin of the much debated soft excess. 

In this paper we have focused on explanations for the 
soft excess which involve blurred reflection from an accre- 
tion disc. Such models have been successful in fitting the 
spectra of previous observations of Mrk 335 (e.g. Ballan- 
tyne et al. 2001; Crummy et al. 2006). In the case of the 
recent Suzaku and XMM-Newton observations, we find that 
two reflectors are required in order to simultaneously explain 
the strong, smooth soft excess and the moderately broad Fe 
line. Our best-fitting model comprises a power law, an in- 
ner, heavily blurred, ionized reflector (which produces most 
of the soft excess) and an outer, almost neutral reflector 
(which is responsible for most of the Fe emission). 

The requirement for more than one reflector is a con- 
sequence of several plausible disc/corona configurations. A 
power-law source in the form of a patchy corona, as well 
as a coronal geometry that changes with time (e.g. due to 
magnetic processes), would for instance give rise to different 
parts of the disc having different ionization states. However, 
about half of the reflected emission in our model originates 
outside about 40 r-g , which seems implausible for a disc with 
a patchy corona. A more centrally concentrated power-law 
source could produce such a situation if it for example con- 
sists of two sources at different heights, so that the inner 
and outer parts of the disc are predominantly irradiated by 
the lower and upper sources, respectively. Alternatively, the 
emission from a power-law source located very close to the 
black hole could be beamed along the disc (e.g. due to rota- 
tional motion about the black hole), which in combination 
with flaring of the disc at large radii would lead to a higher 
illumination of the outer disc. We also note that models 
with multiple reflectors have been successful in modelling 
the spectra of several other Seyfert 1 galaxies, e.g. IH 0707- 
495 (Fabian et al. 2002) and PG 1211-1-143 (Crummy et al. 
2006). It seems likely that such models could be important 
in other sources that, like Mrk 335, exhibit a strong soft 
excess in combination with a moderately broad Fe line. 

It is important to note that our two-reflector model 



for Mrk 335 relies on the fact that the inner reflector with 
^ ~ 300 erg cm produces a very weak Fe line. This is 
because the majority of the Fe in this ionization state is in 
the form of L-shell ions, which in the reflion model are 
assumed to produce very little Ka emission due to resonant 
Auger destruction. However, the amount of resonant Auger 
destruction depends on the conditions in the accretion disc, 
and it is possible that the fluorescent yield for several Fe 
L-shell ions actually remain high (Liedahl 2005). It is not 
clear to what extent this would affect our conclusions. 

It should also be noted that the reflection model used 
for the analysis in this paper assumes a constant density 
throughout the disc. It has been pointed out by Done & 
Nayakshin (2007) that the soft excess is much weaker if a 
hydrostatic model for the disc is used. This is because the 
partially ionized material that produces the soft excess can 
only exist in a thin layer in a hydrostatic disc. However, 
it is not at all clear that accretion discs are in hydrostatic 
equilibrium, as e.g. magnetic fields are likely to be impor- 
tant. This therefore does not pose a serious problem for the 
interpretation that the soft excess is due to refiection. 

Our analysis of the spectral variability of Mrk 335 in 
terms of the two-reflector model shows evidence for a pivot- 
ing power law and suggests that both refiectors are varying 
together with the continuum (although there is some evi- 
dence that the inner refiector is slightly less variable). This 
behaviour is rather different from that observed in the typ- 
ical fiux states of many other refiection-dominated sources, 
such as MCG-6-30-15 (Miniutti et al. 2007), NGC 4051 
(Ponti et al. 2006) and IH 0707 (Fabian et al. 2004). In these 
sources the spectral variability can be described in terms of 
a power law that varies in normalization (but not in slope) 
and a much less variable refiection component. The spectral 
variability in Mrk 335 is also unusual in that the biggest 
changes occur in the low fiux states, where the source hard- 
ens substantially. The reason for this unusual behaviour is 
unclear and more data when the source is in a low state is 
needed in order to properly study this variability. 

As an alternative to the refiection model we have also 
explored the possibility of the soft excess being produced 
as a result of smeared absorption. We have seen that the 
swiNDi absorption model provides a very good fit to the 
spectrum of Mrk 335 as long as we only consider the data up 
to 10 keV. However, the model underpredicts the PIN data 
in the 14-40 keV range, suggesting that more refiection is 
present in the spectrum. The observed fiat rms spectra of the 
XMM-Newton observation and part of the Suzaku observa- 
tion also disagree with the absorption model, as this model 
predicts enhanced variability in the 0.8-2 keV range, assum- 
ing that the ionization state of the absorber responds to the 
changing luminosity of the source (Gierlihski & Done 2006). 
As previously mentioned, the swindi absorption model is 
also problematic in its simplistic treatment of the velocity 
field, and it seems like a more realistic model cannot produce 
the smooth shape of observed soft excesses unless extreme 
velocities are incorporated (Schurch & Done 2007). An ab- 
sorption origin for the soft excess in Mrk 335 hence seems 
very unlikely. 

We have seen that the shape of the hard spectrum of 
Mrk 335 can be very well explained with a blurred reflection 
model. However, it should be pointed out that partial cover- 
ing models have also been seen to provide acceptable fits to 
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the 3-10 keV energy band in e.g. the XMM-Newton obser- 
vations of the source (Longinotti et al. 2007; O'Neill et al. 
2007). We have not considered partial covering models in our 
analysis but we note that the rapid variability of Mrk 335 im- 
plies that the emission region must be very small (less than 
about one light hour across), and a partial coverer would 
hence have to be even smaller and very dense. We also note 
that a broken power law (or some other additional compo- 
nent) is required to explain the soft excess within a partial 
covering model (Grupe et al. 2007), making such a model 
less physically plausible. 

Our comparison of the XMM-Newton and Suzaku ob- 
servations has provided interesting information about the 
origin of the narrow emission lines seen in this object. The 
6.4 keV line is clearly detected in both observations and the 
flux in the line is consistent with being constant, confirming 
that the line originates in distant material, e.g. the molecu- 
lar torus of AGN unification models. The 7 keV line which is 
clearly seen in the XMM-Newton spectrum is, however, not 
detected in the Suzaku data. O'Neill et al. (2007) showed 
that the line cannot be attributed to the blue horn of the 
broad Fe line, and our own spectral and rms analysis support 
this interpretation. Specifically, our rms analysis shows that 
the line does not vary during the XMM-Newton observation, 
which places an inner limit on the origin of the line at about 
a light day from the nucleus. On the other hand, the fact 
that the line has varied in the six months between the two 
observations places an outer limit on its origin at around 
a light month. These constraints together suggest that the 
line originates somewhere in the Broad Line Region of the 
source. 



7 SUMMARY 

In this paper we have presented an analysis of a long 
Suzaku observation of the NLSl galaxy Mrk 335. We have 
also considered a previous XMM-Newton observation of the 
source. 

The Suzaku spectrum exhibits a broad Fe line and a 
strong soft excess, and does not show any signs of warm ab- 
sorption, in agreement with previous observations. We find 
that a model consisting of a power law and two reflectors 
provides the best fit to the entire 0.5-40 keV time-averaged 
spectrum. In this model, a heavily blurred inner refiector 
produces most of the soft excess, while an outer refiector 
(outside ~ 40 rg) produces most of the Fe line emission. 

The spectral variability of the Suzaku observation is 
characterised by spectral hardening at low count rates. Fit- 
ting of different fiux states with the two-refiector -I- power- 
law model suggests that this hardening is mainly due to 
pivoting of the power law. The two refiectors appear to vary 
with the continuum, although there is some evidence for the 
inner refiector being less variable. The rms spectrum of the 
entire observation has the curved shape commonly observed 
in AGN, but the shape is significantly flatter when an inter- 
val which does not contain any deep dips in the lightcurve 
is considered. This suggests that the main driver behind the 
curved shape of the rms is the spectral hardening which oc- 
curs in these dips. The physical reason for this hardening is 
currently unclear. 

The time-averaged XMM-Newton spectrum can be fit- 



ted with a similar two-reflector model as the Suzaku data 
and the two data sets have similar fluxes. However, the 
XMM-Newton data do not show a correlation between the 
spectral hardness and the count rate, and the rms spec- 
trum is flat. Since the spectral hardening of the Suzaku data 
mainly occurs in the deep dips in the lightcurve, it seems 
plausible that the lack of significant spectral variability in 
the XMM-Newton observation is due to the fact that the 
count rate never drops to a low enough level. 

We have discussed the implications of our findings for 
the origin of the soft excess in other Seyfert Is, and we 
suggest that the two-reflector model might be applicable 
in other sources that show a soft excess which is stronger 
and/or smoother than expected from the Fe line properties. 
We have also noted that the observed flat rms spectra, as 
well as the high-energy data from the PIN detector, dis- 
favour an absorption origin for the soft excess in Mrk 335. 
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